Fragility of the glassy state of a directed polymer in random media: 

a replica approach 
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We study the stability of the glassy equilibrium state of a directed polymer in random media 
against weak perturbations including the intriguing problem temperature-chaos. The statistical 
correlation of two real replica systems under perturbations, such as small difference of temperature, 
is studied using the mapping to a quantum problem of interacting bosons. Explicit replica symmetry 
breaking takes place being induced by the infinitesimally weak perturbations in the n — ♦ limit 
followed by the thermodynamic limit which leads to complete decorrelation of the two real replicas. 
The overlap lengths anticipated by phenomenological scaling arguments are retrieved as a correlation 
length associated with this explicit replica symmetry breaking. 



A standing problem of glassy systems with random- 
ness and frustration is the possible instability of the 
glassy frozen states against infinitesimally weak pertur- 
bations such as an infinitesimal change of temperatures. 
A class of phenomenological scaling theories started first 
in the context of spin-glasses [Q suggests that the equi- 
librium glassy state dramatically changes beyond the 
so called overlap length. Unfortunately, the validity of 
the prediction has not been proven explicitly except for 
some Migdal Kadanoff real space renormalization-group 
(MKRG) studies [|. The problem remains very far 
from being resolved especially concerning the problem 
of temperature-chaos, i. e. the sensitivity to small tem- 
perature changes || || . 

The problem of temperature-chaos is now of great in- 
terest since it is considered as a possible mechanism of 
the rejuvenation (restart of aging) found in experiments 
such as spin-glasses || and polymer glasses || . Concern- 
ing the memory effect which is simultaneously observed 
in experiments, recently a novel dynamical memory was 
found in a coarsening model under explicit cycling of tar- 
get equilibrium states [0. 

We consider a much simpler model, the directed poly- 
mer in random media (DPRM) in 1+1 dimension || ||] 
which allows a much better access to the problem. It is 
known to be in the frozen phase at all finite tempera- 
tures in the sense that its scaling properties are always 
governed by the T = glassy fixed point. In spite of its 
simplicity compared to spin-glass models, DPRM is be- 
lieved to possess many of the subtle properties of glassy 
systems. Indeed the anomalous sensitivity of the glassy 
states to weak perturbations are predicted by scaling 
arguments which are partially supported by numerical 
studies Jh], H@||0§. In the present paper, 



we present a unified analytical approach to deal with a 
variety of perturbations based on the concept of explicit 
replica symmetry breaking |l6|, |l^] and confirm analyti- 
cally the predictions of the scaling arguments for the first 
time. 

DPRM belongs to the class of elastic manifolds in ran- 
dom media Jl8| so that its statistical mechanical proper- 
ties are potentially relevant for a broad class of physical 
systems of much interest. The latter includes the domain 
walls of ferromagnets JTs| ] with weak bond randomness 
and the flux lines in type-II superconductors with ran- 
domly distributed point like pinning centers ^0|, CDW 
and vortex lattice systems with weak random-periodic 
pinnings 0, ^2| . 

We study the model described by the following Hamil- 
tonian, 



H[V,h >( f>}= f 1 
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The scalar field cf> represents the displacement of the elas- 
tic object at point z in a 1-dimensional internal space of 
size L. The 1st term in the Hamiltonian is the elastic 
energy, k being the elastic constant. The random pin- 
ning medium is modeled as a quenched random poten- 
tial V{4>, z) with zero mean V((f>, z) = and short-ranged 
spatial correlation, V<j>, z)V(<p', z') = 2DS(<fi~<fi')5(z—z'). 
Here and in the following 777 stands for the average over 
different realizations of the random potential. 

Our starting point is a system of two real replicas, say 
A and B, whose configurations 4>a(z) and </>s(z) are sub- 
jected to exactly the same random potential and temper- 
ature. Second, we apply small perturbations to them. In 
the present letter, we consider 

i) Temperature difference Jig ]: slightly different tem- 
peratures T A = T + ST T B = T - ST for A and B 
respectively with 5T/T < 1 

ii) Decorrelation of random potential JToj , O]: the 
random potential of B is made from that of A 
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as V B = (V A + 6V')/VTTP where \S\ < 1 
and V' follows the same Gaussian distribution as 
V. Then Vg{4>> z) = and V G 4>, z)V G -(<f>', z<) = 
2D GG 'S(<t> - 4>')5 (z - z ') with D AA = D BB = D 
and Dab = D/^/l + 5 2 < D. 

iii) Explicit short-ranged repulsive coupling |l^, ^7): A 
and B replicas are subjected to explicit repulsive 
short-ranged interaction e dz5(4>A{z) — 4>b(z)) 
with < e < 1. 

iv) Tilt field fl^| : A and B replicas are subjected to 
a tilting field of opposite sign —h(j>A(L) + /i</>b(L) 
with h < 1. 

To study the consequence of perturbations, we analyze 
disorder average of partition function Z A+B (L) in which 
both real replicas G = A, B are further replicated into 
n- replicas a — 1, . . . ,n. As noticed by Kardar if 
an analytical continuation for n — * is possible, such a 
replicated partition function can be identified as a genera- 
tor of cumulant correlation functions of sample-to-sample 
fluctuations of free-energies p3] . 



\nZ n A+B {L) = n[-f3 A F A {L)- 13 B F B {L)} 



+ —[-PaFa(L)-PbF b (L)}1 



(2) 



where [. . . ]£ stands for p-th cumulant correlation func- 
tions of the total free-energies —/3aFa(L)—(3bFb(L) with 
Fq{L) and (3 G being free-energy and inverse temperature 
of each subsystem G= A, B. 

Then if there is a complete change of the free-energy 
landscape due to infinitesimally weak perturbations of 
strength { « 1, the statistical correlations between A 
and B should be lost asymptotically in the triple limits 
: n — > followed by L — > oo and finally 5 — ► 0. Here 
the order of the limits is crucial. Then the total partition 
function Z\ +B (8, L) under perturbation should reduce to 
a product of the partition functions of the sub-systems 
Z 7 a(L) and Z B (L) for A and B respectively, 



lim lim 

<5->0L 



hm Z\\ B (S, L) = Z\ (L) x Zl (L) (3) 

oo n— >0 



We use this factorization as the definition of the decor- 
relation of the free-energy landscape (chaos). 

To be specific, we consider that one end of all replicas 
is fixed at x — and the other end fixed at {xG t a}- The 
partition function of the total system can be expressed 
by a path integral over all possible configurations of 2 x n 
replicas, 



Zl +B {{x G , a },L)= Yl V<j} G , a e 



-Sa+b[<I>c 



(4) 



G=A,B 



with a certain effective action SA+B[4>G,a\- Then in turn, 
the partition function can be generated as a path in- 
tegral of a quantum system in imaginary time with a 



Schrodinger equation, 
d 



!L ZJ +B ({x Gta },L)=HZ- +B ({x Gi01 },L). (5) 
with the Schrodinger operator 

" = -£ 



2k dx 2 r , „ 

G,a G > a 

Dec 
(k B T G )(k B TG>) 



E 

(G,a), 
(G',/3) 

+€^S((/)A, a (z) - <f>B,a(z)) 



$(4>G,a(z) ~ 4>G',p{z)) 



E 



d 



h 



k B T dxA, a k B T 



E 



dx 



B.a 



(6) 



Here we have two different kinds of bosons for A and B 
each of which has n particles. The 1st term represents 
the kinetic energy. The 2nd term stands for attractive 
short-ranged interactions between the bosons where the 
sum is taken over all possible pairs of bosons. The 3rd 
term is the explicit repulsive coupling. The last term is 
due to the tilt field applied to A and B subsets in the 
opposite directions. 

The Shcrodinger operator can be decomposed as 



Ti-A+B =7io + SH 



(7) 



where 



^ = -E 

G, a 



k B T d 2 
2k dx 2 - 



D 



G.a 



[k B Tf 



^2 6(xG,a- XG',p) (8) 



((G,a), 
(G',/3)) 



is the Schrodinger operator in the absence of perturba- 
tions and STC is a perturbation term representing one of 
the specific perturbations i)-iv). Note that the original 
operator 7io has a high replica (permutation) symmetry: 
it is symmetric under all possible permutations of bosons 
not only among A and B subsets but among the whole 
set. 

Obviously, the perturbation term 5TL breaks the replica 
(bosonic) symmetry explicitly. Then the interesting 
question is whether such a term can break the symmetry 
in the thermodynamic limit even if it is infinitesimally 
small. If a decorrelation of the free-energy (chaos) takes 
place, it should manifest itself as reduction of this high 
symmetry, i. e. replica symmetry breaking, such that 
replica symmetry remains at most only within each sub- 
set A or B. Our problem can be considered as a particu- 
lar example of the problem of explicit replica symmetry 
breaking proposed by Parisi and Virasoro , who tried 
to give a prescription to define replica symmetry breaking 
in a sound thermodynamic sense by introducing infinites- 
imally weak replica symmetry breaking terms. Indeed 
our approach is an extension of the work by Parisi uM 
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who implemented the idea of explicit replica symmetry 
breaking in the context of DPRM. 

The ground state of the Schrodinger operator || can be 
obtained via the Bethe ansatz ^6|. It is a bound state 
involving all particles described by the wavefunction 



< *Rs|{a:G,a} >~ exp 



51 l x "' G 



X/3,G' 



with A = kD / (feeT) 3 . Since it is completely replica sym- 
metric, we call this as RS (replica symmetric) state. Fol- 
lowing Parisi p7| , we consider also an excited state 



< *rsb| 



=< * A 



RSI 



RSI 



< *Rsl*RS >= 



0. 

(10) 



where < and < ^fgl are similar bound states de- 

scribed by the wavefunction like (j|) but constituted by 
only A and B particles. Assuming that the centers of 
mass of the two bound states are infinitely far from each 
other so that their overlap is zero, such a state becomes 
also an eigenstate. Since replica symmetry is reduced 
here, we refer to the state as RSB (replica symmetry 
broken) state. 

Let us consider the difference (gap) of the eigen value 
of the RS ground state and the RSB excited state under 
the perturbation 



T(n)=T {n)+ST{n) 



(11) 



where To(n) = k^T '/ 'k A 2 n 3 is the original 'energy gap' 
and ST(n) is the correction due to the perturbation. 
Assuming that the amplitudes of the perturbations are 
small enough, one may compute the correction term at 
1st order of perturbation as, 



ST(n) 



;*rsb|(5W|*rsb> (*rs|<SW|*rs) 



(*rsb|^rsb) 



(*RS|^RS) 



(12) 



As we outline later, the correction terms can be com- 
puted for the cases that we consider and we find a generic 
form for the leading term 



-ST(n) ~ An p 1 < p < 3 



(13) 



with A > being a constant which represents the 
strength of perturbations and p being a positive integer 
smaller than 3. We call the exponent p as order of pertur- 
bation. Here 'leading' term means a term with smallest 
power of n which becomes most important in the n — > 
limit that we consider. Let us note here that we have cho- 
sen the detailed form of the perturbations i)-iv) such that 
this gap (13) remains invariant under exchange A <-> B. 
Combining the results one finds, 



r(n) = r (n) 



n* (A) 



-(3-pY 



(14) 



with 



k (A) = 



A 



X 2 k B T/ K 



i/(3-p) 



(15) 



Let us recall that we have to take the n — > limit before 
L — > oo. Since p < 3 holds for all the cases we con- 
sider, an arbitrarily small perturbation A will induce a 
level crossing at n*(A) below which the contribution of 
RSB excited state to the replicated partition function be- 
comes larger than that of the original ground state (RS 
i. e. there is replica symmetry breaking. The result 
strongly suggests that the factorization of the partition 
function (|J), i. c. complete change of the free-energy land- 
scape (or chaos) , takes place in the thermodynamic limit 
even with infinitesimally small (but non-zero) perturba- 
tion of strength A -C 1, 



IS), 
& 



lim lim lim 

A^O L^oo n/n*— >0 



ZT +B (A,L) 



Z\{L) x Z%(L) (p<3). 

(16) 



Now let us further exploit from the above result to 
find a more physical picture. In the absence of perturba- 
tions, the logarithm of the replicated partition function is 



known to have a functional form J24| |25| In Z^ +B (A, L) — 

-0fL(2n) + g(2nL 1 / 3 ) which yields [-/?F]? cx U> e with 
9 being the exact stiffness exponent = 1/3. On the 
other hand, ( |l4"| ) implies n/n* is another natural variable 
of the replicated partition function . 

Combining the two, we conjecture the following scaling 
ansatz, 



]nZf +B (L)+pfL(2n) 



: g(2nL 1/3 ,n/n*) 
:g{2nl}/\L/L*) 



(17) 



where / is the average free-energy per unit length. Here 
we introduced a characteristic length L* defined as 



L* 



(71*)- 



A -3/(3-p) 



(18) 



The n — > limit induces the thermodynamic limit 
L — > oo if the variable nL 1 ^ 3 = x is fixed. Then for 
a fixed x, two limiting behaviors are expected. First 
the small length scale regime L/L* <C 1 corresponds 
to n/n* > 1 so that the effect of perturbation will 
be small and the partition function is essentially the 
same as that of the unperturbed system of 2 x n-replicas 
g (x,L/L* — > 0) ~ g(2x). Second the large length regime 
L/L* ^> 1 corresponds to n/n* <C 1. Thus the decou- 
pling ( p^| ) implies g(x, L/L* — > oo) ~ 2 x g(x). 

The above result implies the crossover length L* should 
correspond to the overlap length L c anticipated by the 
scaling argument. In the following we present the outline 
of our results for the specific perturbations. Details will 
be presented elsewhere l2l|. 
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Temperature and potential change 

In the case of i) small temperature difference, one finds 
by the perturbation calculation that the correction to the 
gap ( |l3| ) is characterized by order p = 2 with amplitude 
A oc {5T/T) 2 . Then using (JjJ) and @, we find the 
overlap length L* ~ (8T/T)~\ Similarly, the case of ii) 
small decorrelation of random potential also amounts to 
the same order p — 2 with A ~ <5 2 . Then one finds the 
overlap length L* ~ <5~ 6 . Quite remarkably, the overlap- 
length obtained here is in agreement with the prediction 
of the scaling arguments |0, p"5) . Furthermore, it is in- 
teresting to note that the two perturbations which look 
very different at a first sight, turn out to be the same in 
the replica space. 

Explicit repulsive coupling 

The case of hi) explicit short-ranged repulsive coupling 
has been considered by Parisi [J17J and the result gives 
order p = 1 and amplitude A ~ e. Thus one finds the 
overlap-length L* ~ e~ 3 / 2 which agrees with the predic- 
tion Q based on a scaling argument. 

Tilt field 

Finally we consider the case of iv) tilt field. In this case, 
the 1st order correction term (13) vanishes. Fortunately, 
we can proceed as the following. The Schrodinger oper- 
ator with the tilt field can be rewritten into the fully 
symmetric form as the unperturbed one (|J) by shift- 
ing the momenta, d/dx' Aa — d/dxA.a — h/k^T and 
d/dx' Ba = d/dxB, a + h/ksT. Therefore the exact 
ground state under the tilt field can be obtained from 
the Bethe's wave-function (^) by undoing the previous 
shifting of moments. On the other hand, the unper- 
turbed RSB wavefunction ( |l0| ) remains as a valid eigen- 
state under the tilt field. Combining these results, we 
obtain the correction to the gap due to the tilt field as 
ST(n) = nti 1 / nk^T from which we read off order p = 1 
and amplitude A ~ h 2 , which yields the overlap length 
L* ~ hr z . The latter agrees with the prediction |Q 
based on a scaling argument. 

To summarize we have presented a unified analytical ap- 
proach to the problem of the sensitivity to perturba- 
tions of the equilibrium glassy states of 1+1 dimensional 
DPRM based on the idea of explicit replica symmetry 
breaking. Our result in replica space suggests that per- 
turbations can be classified according to their order p and 
the symmetries that are preserved. Therefore, apparently 
different perturbations can bring about the same effects. 
The universal features of the crossover from the weakly 
perturbed regime L/L* < 1 to the strongly perturbed 
regime L/L* 1 can be examined numerically |2j|. The 
prediction of the scaling arguments is supported for 4 
different kinds of perturbations. It may be promising to 
apply the approach of the present paper in the studies of 



other glassy systems including spin-glasses. 
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